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Abstract. Although there is overwhelming evidence of dark matter from its gravitational 
interaction, we still do not know its precise gravitational interaction strength or whether it 
obeys the equivalence principle. Using the latest available cosmological data and working 
within the framework of ACDM, we first update the measurement of the Newton’s constant 
for all matter: X 10 m^kg ^s which differs by 2.2a from the standard 

laboratory-based value. In general relativity, dark matter equivalence principle breaking can 
be mimicked by a long-range dark matter force mediated by an ultra light scalar field. Using 
the Planck three year data, we find that the dark matter “fifth-force” strength is constrained 
to be weaker than 10“^ of the gravitational force. We also introduce a phenomenological, 
post-Newtonian two-fluid description to explicitly break the equivalence principle by intro¬ 
ducing a difference between dark matter inertial and gravitational masses. Depending on the 
decoupling time of the dark matter and ordinary matter fluids, the ratio of the dark matter 
gravitational mass to inertial mass is constrained to be unity at the 10“® level. 
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1 Introduction 

Although Newton’s description of gravity has existed for over three hundred years, the value 
of the Newtonian gravitational constant Gm has only been measured at the 10“^ level. This 
is in contrast to the hne-structure constant, which has been measured very accurately with 
a relative error of 10“®. Over the past thirty years, there has been a large effort to precisely 
measure the value of the Newtonian gravitational constant in the laboratory. Most of these 
measurements were torsion balance based and achieved a level of accuracy which can constrain 
many new physics scenarios (see Ref. [1] for a review). However, these laboratory based 
measurements suffer from systematic errors and have results which are inconsistent with 
each other. As a result, the discrepancy between the measurements lead the Committee on 
Data for Science and Technology (CODATA), which determines the internationally accepted 
standard values, to increase the relative uncertainty of the 2010 recommended value of Gv 
by 20% to 120 parts per million compared to that of the 2006 value [2]. Additionally, a recent 
measurement using cold atom interferometry obtains similar accuracy [3], but the central 
value is in tension with the CODATA standard value at 1.5a. 

In addition to laboratory-based measurements of Gv, cosmological data can also provide 
an independent measurement at a much larger distance scale and also serve as a consistency 
check. As pointed out in Ref. [4] and employed in Refs. [5, 6], the measurements of the 
Cosmic Microwave Background (CMB) temperature anisotropy and polarization can be used 
to extract the value of Gv from cosmology. In the last few years, an enormous amount of 
cosmological data both from satellite and ground-based experiments has been gathered. The 
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Planck satellite in particular has ushered in the era of precision cosmology, measuring the 
CMB temperature and polarization anisotropy with unparalleled accuracy [7, 8]. In addition, 
ground-based telescopes such as the Atacama Cosmology Telescope (ACT) [9] and the South 
Pole Telescope (SPT) [10], provide high accuracy polarization measurements at high multipole 
moments. In this work, we will update the cosmological measurement of Gn for all matter 
using the latest available cosmological data, including the 2013 Planck release, ACT/SPT, 
and Big Bang Nucleosynthesis (BBN). 

Another interesting question that may be answered using cosmological data is whether 
dark matter has the same gravitational interaction as the ordinary baryonic matter. So far we 
have seen extensive evidence for the existence of dark matter, such as galactic rotation curves 
and gravitational lensing. There is no doubt that dark matter has gravitational interactions, 
although we have not found additional forces felt by dark matter, despite a large effort to 
understand dark matter particle properties. This serves as motivation to understand the dark 
matter gravitational interaction better. In this paper, we will use the current cosmological 
data to constrain the properties of the dark matter gravitational interaction. 

Different gravitational interactions for different matter types inevitably break the Weak 
Equivalence Principle (WEP), which states that all objects in a uniform gravitational field, 
independent of the mass or other compositional properties, will experience the same acceler¬ 
ation. In Newtonian language, the difference between inertial mass (the mass appearing in 
Newton’s second law) and gravitational mass (the mass appearing in Newton’s law of gravity) 
must be exactly zero for the WEP to be respected. For ordinary baryonic matter, modern 
experiments using torsion balances report that the difference between inertial and gravita¬ 
tional masses is zero at the 10“^^ level [11]. Thus, violations of the WEP in the visible sector 
are tightly constrained. The constraint on WEP breaking in the dark sector is much less 
restrictive. 

One simple way to mimic dark matter equivalence principle breaking is to introduce 
a long-range force only for dark matter [12-16]. In this class of models, the gravitational 
interaction strengths for dark matter and baryonic matter are the same, although the total 
long range force between two clumps of dark matter is different from baryonic matter. In our 
study, we will also update the cosmological constraints on the fifth force strength and pay 
additional attention to the force carrier energy density. 

Beyond the fifth-force description, we also introduce a phenomenological, post-Newtonian 
model to explicity distinguish dark matter gravitational and inertial mass. For ordinary cos¬ 
mology, one can derive parts of the background Friedmann and linear perturbation equations 
in the post-Newtonian description by treating the scale factor as an expanding sphere filled 
with a homogeneous and isotropic fluid [17]. For pressureless matter, the continuity and Euler 
equations of fluid mechanics can be matched to the energy-momentum conservation equations 
in general relativity. We will follow this description and treat baryonic matter and dark mat¬ 
ter as two separate but coupled fluids. Effectively, the dark matter fluid is evolving with a 
different scale factor such that the ratio of dark matter gravitational to inertial mass can not 
be simply absorbed by the dark matter energy density. The cosmological data will provide a 
constraint for this ratio. A relevant example from the literature uses the tidal disruption of 
the Sagittarius dwarf galaxy orbiting the Milky Way to constrain this ratio at the 10% level 
[18]. We will find that cosmological data can provide a much more stringent constraint. 

Our paper is organized as follows. We first constrain Newton’s constant for all matter 
in section 2. Then, in section 3 we update the constraints on dark matter fifth forces and in 
section 4 we develop a two-fluid description to constrain dark matter WEP breaking. Finally 
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we conclude our paper in section 5. 


2 Newton’s Constant for all Matter 

It is well-known that the gravitational or Newton force of a probing body only depends on the 
product of the Newton’s Constant G and the central body mass. To break this degeneracy, 
an additional force (such as the electromagnetic or the weak interaction) is required to define 
the central body mass. Because of this fact, the existing studies in the literature [4-6, 19] 
have used data from the primordial abundances of light elements synthesized by BBN and 
CMB anisotropies to constrain G, or equivalently other fundamental constants like the fine- 
structure constant [20]. In this section, we use the currently available cosmological data to 
derive a constraint on Newton’s constant. 

Before we go into our detailed analysis, we introduce a dimensionless parameter Xq to 
quantify the potential deviation of Newton’s constant G from the value, Gjsf, measured in the 
laboratory-based experiments 


G — Xq Gn . (2-1) 

We will use the currently suggested central value of Gn = 6.67384(80) x 10“^^ m^kg“^s“^ 
from CODATA 2010 [2]. ^ The relative error is 1.2x 10““^. The latest measurement using laser- 
cooled atoms and quantum interferometry reaches the similar precision and has an agreement 
with the CODATA value at 1.5a [3]. The constraints from cosmological data will serve as an 
independent, time-sensitive measurement at large length scales. 

2.1 Dependence of CMB Anisotropy on Newton’s Gravitational Constant 

With the introduction of Xq in Eq. (2.1), the Friedmann equation is modified as 

=(;()' = ( 2 . 2 ) 

where 7i is the Hubble rate of expansion; a is the scale factor of the universe; p is the total 
energy density; dot indicates derivative with respect to the conformal time r with dt = a{T)dT. 
From Eq. (2.2) we see that the effect of Xg is to change the amplitude of the expansion rate 
'H{XG,a) by a factor of Xg, but not the shape. Since no new preferred cosmological scale is 
introduced by varying G, a variation in Newton’s gravitational constant is equivalent to a 
simple rescaling of the wave numbers, as pointed out in Ref. [4]. Using the baryonic equations 
of motion in the conformal Newtonian gauge [21], one has 

Sb = -9b + 3^, (2.3) 

9b = --9b + (?s Sb -h aneaT{9.y - 9b) + ■ (2.4) 

a 3pb 

Here, 6b = Spbjpb', 9b = *^^6’ baryon speed of sound; rie is the free electron number 

density; ax = 0.6652 x 10“^'^ cm^ is the Thomson scattering cross section; (j) and ^ are the 
scalar metric perturbation. If ax was zero or there was no Coulomb interaction, one can 
show that Eqns. (2.2)(2.3)(2.4) are independent of Xg by a replacement of r —>• tXg and 

’There are large inconsistencies among different measurements. The CODATA Task group has taken the 
11 values after each of their uncertainties multiplied by an ad-hoc factor of 14. 
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k —)• k/Xc- The density fluctuations produced by a mode of wave-vector fc in a universe with 
Xg 7 ^ 1 have dynamics equivalent to a mode with fe' = /c/Ag in a universe with Xq = 1- 
Thus, for primordial fluctuations governed by a power law, changing G can be compensated 
by adjusting the amplitude of the primordial power spectrum appropriately. As a result, 
the CMB anisotropy spectrum does not change when varying Xq if there is no Coulomb 
interaction. This makes sense physically because the effect of changing G is to cause the 
universe to expand slightly faster or slower by a factor of Xq, causing the “expansion clock” to 
run at a different rate. Since we only measure angles in the CMB through ratios of distances, 
such a change precisely cancels. 

To observe the effect of different values of G in the CMB, we need an independent 
clock that measures the expansion rate. This independent clock will come from the physics 
of recombination through an additional Coulomb force acting on the baryons. The number 
density, Ug, of free electrons will then depend on Xg- As we will show, this change in Ug is 
observable in the CMB through its effect on the visibility function that enters into the integral 
solution for the temperature anisotropies produced by a mode of wave-vector k observed 
towards direction h [22]. The temperature at a direction h is an integration of the sources in 
the line of sight convoluted by the visibility function (/(r), which is related to the free electron 
number density ng{T) through the Thomson scattering as 

d r° 

5 f(r) = — — [exp(— k)] , with k{t) = ar / a(r)ne(r)dr. (2-5) 

dr 


Here, tq is the current universe time. The CMB anisotropy dependence on the visibility 
function is precisely what makes a change in the gravitational constant detectable. This is 
because the free electron density Ug depends on the physics of recombination through the 
ionization fraction Xg = ngjnn which evolves according to [21] 


dxg 

dr 


a Cr 


m)ii-xg) 


nHa^'^\Tb)xl , 


( 2 . 6 ) 


where Cr is the Peebles correction factor to account the presence of non-thermal Lyman-a 
resonance photons; Tb is the baryon temperature; nn is the number density of hydrogen 
atoms. Here, l3{Tb) is the collisional ionization rate from the ground state 



mgksTb 

2Trh? 


3/2 




(Tb), 


(2.7) 


with Bi = 13.6 eV and a^‘^\Tb) as the recombination rate to excited states 




647r (ksTb 


- 1/2 

MTb), 


02(Tfe) ~ 0.448 In 



( 2 . 8 ) 


Changing variables from r to the scale factor a, d/dr = oHd/da = Xg f{a) ad/da, we 
have Eq. (2.6) in a different form 


dxg 

da 


Cr 

AG/(a) 


P(Tb){l - Xg) - nHa^‘^\Tb)xl 


(2.9) 


The above equation shows that the form of Xg{a) depends on Xg- For a larger value of Xg-, 
this form shows a larger value of Xg for a given a or redshift. This can be understood as that 
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Figure 1. Left panel: the ionization fraction Xe as a function of z for different Ac- Right panel: 
the visibility function ^(t) as a function of conformal time. The axes have been rescaled using the 
recombination time rjj to remove the overall scaling of r with Ac- 


for a larger Ac, the universe expands faster at a given redshift, since it becomes more difficult 
for hydrogen to recombine and this leads to a larger value for Xg- This effect is demonstrated 
in the left panel of Fig. 1 in terms of the redshift z = 1/a — 1. 

Since the ionization fraction determines the density of free electrons, we also expect 
the visibility function to change when varying Aq. The normalized visibility function can 
be thought as the probability that a photon last scattered at a particular conformal time 
r. As Ag increases, recombination (and the last scattering of photons) takes place over a 
longer period of time which means the visibility function becomes broader. This broadening, 
shown in the right panel of Fig. 1, leads to a damping of anisotropies on small scales on which 
photons are still scattering. In Fig. 2, we show the net effect of changing G on damping (for 
Ag > 1) or enhancement (for Ag < 1) of the temperature anisotropies, with an emphasis 
on the small angular scales. Also in Fig. 2, we show the EE and TE power spectra. Since 
the effects are more dramatic at a small scale or a large £, we will later use this fact to 
understand constraints from different experiment data. We also note that the effects from 
varying Ag have a large correlation with the parameter Ug of the primordial power spectrum, 
since introducing an appropriate tilt of this spectrum can also act to damp or enhance the 
small scale peaks. 

2.2 Analysis Method 

We perform a Markov Chain Monte Carlo analysis using the publicly available Monte Python 
code [23] which interfaces with the CLASS code [24, 25]. In addition to Ag, we sample 
the concordance ACDM parameters which include the physical baryon and CDM densities, 
Ub = and iOc = the Hubble parameter Hq at the current time, the scalar spectral 

index n^, the primordial power spectrum normalization In(lO^^As) at fco = 0.05/Mpc and the 
reionization optical depth r. Flat priors were used on all the above cosmological parameters. 

The chains are checked for convergence using the Gelman-Rubin R—1 statistic for each 
parameter [26]. To obtain constraints on the cosmological parameters, the Monte Python 
package marginalizes over the remaining nuisance parameters. For computing the likelihood 
we use the package provided by the Planck team with the 2013 data release [27]. It contains 
high and low-£ TT likelihoods in addition to low-.^ TE and EE likelihoods from WMAP9. 
Also included are high-^ TT likelihoods using 3 years of data from the Atacama Cosmology 
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Figure 2. The effects of varying Xq on the TT, EE, and TE power spectra. 


Telescope (ACT) [9], the South Pole Telescope (SPT) [10] and a Planck leasing likelihood. 
ACT measures the CMB angular power spectra over a 600 square degree patch of sky at 148 
and 218 GHz. SPT does the same measurement over a 800 square degree patch of sky at 
95, 150, and 220 GHz. The combined ACT/SPT package covers a multipole moment range 
500 < i < 3500 for use in constraining cosmological parameters. 

In addition to the Planck package, we use likelihoods for the Baryon Acoustic Oscillations 
(BAO) and the Hubble Space Telescope (HST) available with the Monte Python code. The 
HST likelihood comes from [28], which determines the Hubble constant using the Wide Field 
Camera 3 on the Hubble Space Telescope to observe over 600 Cepheid variables in the host 
galaxies of 8 recent Type la supernovae in the optical and infrared. The dataset covers a 
redshift range of 0.01 < z < 0.1. The BAO package contains data from the Sloan Digital 
Sky Survey (SDSS) (Data Releases 7 and 9) [29, 30] and the Six degree Field Galaxy Survey 
(6dFGS) [31]. These experiments have a mean redshift of z ~ 0.10 {z ~ 0.5) for SDSS Release 
7 (9) and 2 : ~ 0.05 for 6dFGS. SDSS Release 7(9) covers 11,663(14,555) square degrees and 
6dFGS covers ~ 17,000 square degrees of sky. 

2.3 Constraints on G 

Using the dataset Planck+ACT/SPT+Lensing+BAO+HST, we report in Table 1 the con¬ 
straints obtained on the cosmological parameters which were sampled using Monte Python. 
This dataset provides a constraint on Xq at about the ~ 2.2% level. Including the BBN data 
raises the mean value of Xq and tightens the constraint to the ~ 1.8% level. To show correla¬ 
tions between Xq and some common cosmological parameters, we show the one and two sigma 
confidence contour plots on the Xq — ng planes in Fig. 3. The degeneracy with the scalar 


6 
















































spectral index is to be expected since varying Ac either delays or hastens recombination, 
which damps or enhances small angular scale oscillations. 
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Table 1. Constraints on the 7 free cosmological parameters (ACDM + Aq) nsing the experimental 
dataset Planck+ACT/SPT+Lensing+BAO+HST. 
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Figure 3. 68% and 95% likelihood contonr plots on the ng—Xc plane using the 

Planck+ACT/SPT+Lensing+BAO+HST dataset. 



We show the constraints on Xq from different combinations of datasets in Table 2. 
Comparing constraints from different groups, one can see that they are consistent among 
each other. All of them have the central values to be above unity at a slightly over one sigma 
conhdence level. The datasets of Planck+ACT/SPT with more data at high-^ provides a 
strong constraint with the relative error of ~ 2.4% on Xq- So, we report the cosmological 
measurement of Newton’s gravitational constant using the combination of Planck, ACT/SPT, 
Lensing, BAO, HST and BBN as 

G (cosmological) = A^ GAr(CODATA) = 7.261q'27 x 10“^^ m^kg“^s“^ , (2.10) 

which has a relative error of 3.7%. This cosmologically measured value is roughly consistent 
with the CODATA and has a 2.2cj tension. 
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Table 2. Cosmological measurement of Xc (with one-sigma errors) from different combinations of 
data. 


3 Scalar-mediated Long-range Dark Matter Force 


The conservative way to differentiate the long-range forces among dark matter and ordinary 
matter is to introduce an ultra light scalar mediator, which only couples to dark matter 
particles. There is existing literature on constraining additional Yukawa-like interactions for 
dark matter [12-16]. Following the notation in Ref. [12], the interaction Lagrangian for the 
dark matter and the scalar mediator is 

£ D ^xXX- (3.1) 

Here, the ultra-light pseudo Nambu-Goldstone boson (j) could be associated with spontaneous 
breaking of some global symmetry at the scale dvr/ [32]. In principle, more complicated terms 
could enter the / potential F(<^). For our later purpose, the (p field value is small such that 
only the leading mass term in the potential is important. For simplicity, we assume that (p 
does not couple to the ordinary matter. 

Mediated by the new light scalar field, two dark matter particles with masses rriDi 
171^2 have the following static potential 

V{r) = _ ^NmD,mD 2 ^ ^ (3 2 ) 

where df = f‘^ and Mpi = I/^/SttGn « 2.4 x 10^® GeV is the reduced Planck scale. 

The light scalar field can also contribute to the energy-momentum tensor and therefore 
can modify the Friedmann equation 






(j) I 


+ V{P) 


(3.3) 


Here, p is the total energy density and pc is the cold dark matter energy density. The 
homogenous scalar p = />(r) has the following equation of motion 


p + 2% p + 
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3.1 Conditions to Ignore the cj) Background 

To simplify our discussion for constraining the dark matter fifth force, we will work in the 
parameter space where we can ignore the (p contribution to the background evolution (ie., 
the Friedmann equation). For an ultra-light light scalar satisfying 


2 1 




(3.5) 


The solution to the equation of motion in Eq. (3.4) is 


where ci is an integration constant and Hq and are the Hubble parameter and dark matter 
density at the current time, respectively. Starting from a radiation-dominated universe with 
= Hq a“^, we have the following solutions (for a particular choice of initial conditions) 


pR 

pR 


l{^) 


" --HoMpi . 


(3.7) 

(3.8) 


In the matter dominated region with 


Hq a ^, we hnd the following solutions 


4>m = —2 



4>m = -2 


^7pi 

/ 




Mpi 


ln(a) 


2C^ 




HqMpi 




(3.9) 

(3.10) 


The two integration constants and can be determined by requiring both p and p to 
be continuous functions: 4>R{aeq) = 4>Mio,eq) and 4>R{aeq) = (/>M(oeq)- The equality of matter 
and radiation happens at Oeq ~ 1/3600, which determines the integration constants 


Cf 



3/2 

eq 


-1.16 X 10"®, 




7 

-InOen ~ 8.77. 

12 ^ 


(3.11) 


To ignore the p background contributions to the Hubble parameter in Eq. (3.3), we 
need to have all three (/>-related terms to be suppressed. It turns out that the requirement, 
p Pel f < /O; in the matter-dominated era provides the most stringent bound and requires 


/ > Y'Tc^Mpi « 4.2Mpi, or a/<0.06. 
The bound on the scalar mass from Eq. (3.5) can be simplihed to be 


(3.12) 


Ho 


< 


cr 


0 . 2 , 


(3.13) 


or rrifj) < 0(10 eV). In the following numerical analysis, we will stay in the parameter 
space satisfying Eq. (3.12) and Eq. (3.13). 
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3.2 Linear Perturbation Equations 

To derive the linear perturbation equations, we expand the scalar held into 


= 4>o{t) + (p{x,T), (3.14) 

with (Po{t) as the background held and (p{x,T) as the hrst order perturbation function. Per¬ 
turbing the (p equation of motion, we arrive at the following equation for the ip evolution [13] 

ip + 2'Hp + {k^ + a^vn?^)p = -]:h(pQ - ^ Pc^cO^ ■ (3.15) 

Here, the function h{k, r) is the scalar metric perturbation in the synchronous gauge [21]. The 
new second-order differential equation in Eq. (3.15) requires two initial conditions for the p 
held. Deep in the radiation epoch, keeping the leading terms in powers of a and neglecting 
k"^ and terms, we have {d/dT){a^ip) = —^ hcpQa?. Using the scaling solutions of /i oc 
and a (X T and noticing that pQ is a constant from Eq. (3.8), we have the following initial 
conditions 

1 ■ 1 - 

p = -—hpoT, p = --hpo. (3.16) 

The hrst-order perturbed Einstein equations are 

k'^ rj — -T-L h = AttGn oP , (3.17) 

k‘^p = AttGn g? ikHT\ , (3.18) 

h + 2nh-2k^r] = -SnGNa^ (3.19) 

h + Qp + 2-H{h + Qp)-2k^p = 2 A-kGn a^ikk - - 5\T\/2,) . (3.20) 

The new contributions to the energy momentum tensor from the scalar held are calculated 
to be 

+ Podc] - - mlpQip, (3.21) 

ik^5T^iiP) = (3.22) 

^T\{p) = ^Pop - SmlpQip, (3.23) 

which agree with the formulas in Ref. [13]. The dark matter density perturbation still obeys 
6c = —^h, simply from the conservation of energy-momentum. Since the held po and p 
change the values of h through the Einstein equations, the dark matter density perturbation 
is affected by pQ in this indirect way. 

3.3 Constraints on af = 

Because of our choice of initial conditions for the scalar mediator evolution, the main con¬ 
straints on the strength of the hfth force come from the perturbation part rather than the 
background evolution. The CMB temperature anisotropy turns out to provide the most strin¬ 
gent constraint. As already mentioned in Ref. [33, 34], the additional long-range attractive 
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force among dark matter particles can introduce the late time variation of the gravitational 
potential. The ISW effect acts to increase the power for low multipoles. To conhrm this 
observation, we show the conformal time-derivative of the gravitational potential in the left 
panel of Fig. 4. Since the universe is expanding, the magnitude of the potential is decaying in 
time and as such the derivative is negative. Comparing to the ordinary cosmology with a/ = 0 
and from Fig. 4, one can see that increasing the value oi af leads to a larger gravitational 
potential and a more dramatic ISW effect. 


xlO-^ 



Multipole Moment ( 



Figure 4. Left panel: the derivative of the gauge-invariant scalar metric perturbation as a function 
of multipole moments, i, for different values of the fifth force strength af and at the current time. 
Right panel: the CMB TT power spectra in i for different values of o/ and m^. The curves for 
= 0 Etnd mcjj/HQ = 0.5 are overlapping each other. Also shown here is the Planck three year 
data points for i < 50 [27] . 


In the right panel of Fig. 4, we show the CMB temperature anisotropies in terms of the 
multipole moments. Since the ISW effect mainly affects the small i region, increasing the 
value of a/ leads to a larger power spectrum. Within the validity of our approximation of 
ignoring the scalar held contribution to the background evolution, the difference between a 
massless m(j,/Ho = 0 and a small mass such as m^/Ho = 0.5 is negligible. So, we do not 
anticipate that the CMB temperature anisotropies can constrain in the regime where the 
approximation condition in Eq. (3.13) is satished. 



Figure 5. The constraints on the fifth force strength aj and the mediator mass from the 
Planck+Lensing data. 


- 11 - 








































































In Fig. 5 and using the data from Planck+Lensing, we show the cosmological constraints 
on the fifth force strength and the mediator mass. One can see that the data is insensitive to 
the mediator mass, but does provide a stringent bound on the fifth force strength as 

a/= < 1.62 X 10“^ (95%C.L.). (3.24) 

Compared to the results in Ref. [16], which used the WMAP 5-year data [35], our results using 
the Planck data show a dramatic improvement on constraining an additional long-range force 
among dark matter. 


4 Dark Matter Weak Equivalence Principle: Two-Fluid Description 


In this section, we use cosmological data to constrain the difference between inertial and 
gravitational mass for the dark matter. This immediately requires a violation of WEP. For 
baryonic matter, modern torsion-balance experiments report that the difference between iner¬ 
tial and gravitational masses is zero at the 10“^^ level [11] (see also [1, 36]). Thus, violations 
of the WEP in the visible sector are tightly constrained. Therefore, in this section we consider 
a violation of the WEP in the dark matter sector only. We introduce a new parameter Xu 
such that the ratio of dark matter gravitational to inertial mass is given by /mD 

with as the inertial mass. For baryonic particles, the gravitational mass is assumed to 
be the same as the inertial mass. For two baryonic particles bi and 62 and two dark matter 
particles Di and D 2 , the gravitational forces in terms of the particle inertial masses are 

GNmh^nib^ „ . GNmb^mDj 2 GNmoimOi .n 

2 ) -^bi,Di ■^D 2 ’ ^Di,D 2 •^D 2 ■ 


E)i,62 


There is a simple relation, ~ -^ 61 , 01 -^ 62 ,£> 2 ) which does not hold in section 3 where 

there is no “fifth force” between one dark matter and one baryonic particles. 

As pointed out in the textbook [17], one can derive parts of the Friedmann and linear 
perturbation equations of cosmology using a post-Newtonian description. The Friedmann 
equation can be derived by considering a sphere filled with a homogenous and isotropic mat¬ 
ter density distribution and studying the radius change as a function of time. For pressureless 
matter, one can use the continuity and Euler equations of fluid mechanics to match to the 
energy-momentum conservation equations from general relativity. In our case with differ¬ 
ent Newtonian forces for dark and ordinary matter, we use this post-Newtonian language 
to motivate modifications to the cosmological background and linear perturbation equations. 
Following the same approach as Ref. [17] takes with a single scale factor for all matter, we 
derive two coupled “Friedmann equations” which govern the expansion of dark and ordinary 
matter as two coupled fluids, each with their own scale factor. We would like to note that we 
are using this post-Newtonian, non-relativistic description in order to motivate a phenomeno¬ 
logical parameterization of dark matter WEP breaking, since true WEP breaking is unlikely 
to be realized within the framework of general relativity. A full, generally covariant model 
that may encode these deviations is beyond the scope of this work. 

The classical picture for our two-fluid description is two expanding spheres with the 
centers located at the same point. The mass densities, pb and pu, for both fluids are assumed 
to be uniform inside the spheres. For a probing baryonic matter particle with mass nib on 
the surface of a sphere with radius r^, its motion is governed by the amount of baryon and 
dark matter inside the Vb radius sphere 

(fn 


nib 


dt'^ 


GNmbMbin) - GNTUbMoin) 
- 9 - “ '^D - 9 - 


(4.2) 
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where Mh{r) = Airr'^pijjS and Mnir) = Here, to describe the particle motion, all 

masses mf,, Mb and Mjy are inertial masses. Similarly, for a probing dark matter particle with 
mass rriD and at a radius r^, one has the equation of motion to be 


ruD 


(frp 

dt"^ 


—Ad 


Gn mp Mb{rp) 


r 


2 

D 


a 2 


Gn mp Mp{rp) 


D 


D 


(4.3) 


In general, the two radius functions in time rb{t) and rp{t) could be independent of each 
other. So, to match to the Friedmann equations, we need to introduce two scale factors, a{t) 
and ap{t), for baryonic matter and dark matter, respectively. 

Requiring the radii proportional to the scale factors, we have rb{t) = r^a{t)/aP and 
= f'^pO-oit)/oPp. We then rewrite Eqs. (4.2)(4.3) as 


1 d?a 
a dt‘^ 


47rG^ 

3 


[pb{a) + Pd PD{aD)\ , 


1 d^ap 
ap dt'^ 


47rG7V 

3 


[Pd Pb{a) + Pp pp{ap)] . 

(4.4) 


Here, the densities are diluted with the expansion of the sphere and are pb = 
and Pp = p^j^[aPj^/ap{tyf. The equations in Eq. (4.4) describe the acceleration of baryonic 
and dark matter particles, when the dark matter sector violates the WEP. The power of Ad 
corresponds to the number of dark matter masses present in the interaction, e.g. a dark matter 
particle moving in response to a baryonic source receives one power of Ad- We see that in the 
limit where the WEP is restored (Ad —^ 1), both equations reduce to the usual Friedmann 
equation as expected. Extending the above analysis to include the radiation energy density 
H/j and dark energy density Ha by assuming they couple to gravity in the same way as the 
baryons, we have our final Friedmann equations 


1 d?a 
a dt"^ 

1 d^ap 
ap dt"^ 


Hi 


Hi 


a3 

Ad 


+ ^ + (1 + 3w^) 


Ha , Pd^d 


q3ui+3 


+ 




Hfe 


+ ^ + (1 + 3-^) 


Ha 

q3w+3 


+ 


A^Hd 


'■D 


(4.5) 

(4.6) 


Here, we introduce the parameter Hq to draw a distinction from the ordinary Hubble constant. 
Ho = ^^Itoday; at the current universe. For \p = 1 with a single scale factor, we have 
Hq = Hq. For the ordinary ACDM cosmology with a flat universe, the dark energy density 
is not an independent parameter and is given by Ha = 1 — H;, — H/j — Hd. For our case with 
two scale factors, we will keep Ha as a free parameter. 

Compared to ordinary cosmology with a single scale factor, we have one additional 
second-order differential equation, which requires two more initial conditions. A simple way 
to hx the dark matter scale factor initial conditions is to have ap = a and da/dt = dap/dt 
at some time a = ajnit. We assume that the dark WEP breaking is turned on at this time, 
parametrized by zp = — 1. Physically, this transition redshift zp corresponds to a scale 

where some interaction coupling dark matter to the Standard Model falls out of equilibrium. 
Before the transition redshift zp, everything evolves as a one-component fluid described by 
a single scale factor a. After zp, the dark matter decouples from the other components and 
evolves as a separate fluid according to a dark scale factor ap. 

Using the measured cosmological parameters from Planck [7], we show the behaviors of 
the two scale factors as a function of time in Fig. 6. In this plot, we choose Oinit = 10“® 
and use the equations of ordinary cosmology to calculate da/dt at the time corresponding to 
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Figure 6. The ordinary matter and dark matter scale factors as a function of time. The cyan and 
dotted line is for ordinary cosmology with Xjj = 1. The two green five-stars correspond to the time 
with the current observed Hubble constant. 

“init = 10 “®. We then treat this as the initial condition to solve the conpled eqnations in 
Eqs. (4.5)(4.6) for two different Xd = 1 — 0.1 and 1 -|- 10“^. As can be seen from Fig. 6 and 
for the A/r < 1 case (red and black lines), both scale factors increase faster than the standard 
cosmology with A^) = 1. The dark matter scale factor increases faster than the ordinary 
matter one. This is dne to the fact that less matter leads to an open nniverse and the 
effective matter for the qd eqnation of motion is smaller than the one for a. For the A^) > 1 
case, a closed nniverse will be obtained even before one obtains the cnrrent measnred Hnbble 
constant. Therefore, we generically anticipate more stringent constraints for the A/r > 1 case 
than the A^i < 1 case. For a larger valne of Oinit or a smaller valne of zt, the difference 

between the new scale factor and the standard one become smaller. This is simply becanse 

the dark matter flnid has less time to evolve deconpled from the baryonic matter flnid. 

4.1 Dependence of CMB Anisotropy on Xd 

So far we have only motivated modihcations to the backgronnd evolntion eqnations. To 
motivate modihcations to the linear pertnrbation eqnations for the dark matter hnid, we will 
work in the baryon co-moving frame and dehne the conformal time nsing the ordinary baryon 
scale factor dt = a{T)dT. The linear pertnrbation eqnations for ordinary baryons stay the 
same. The linear pertnrbations for the dark matter hnid receive the following modihcations 
(see Appendix A for a detailed explanation) 

+ {Td — T-Id) {^ + kdk) 6 d = 3 6(j} — 9D, (4.7) 

9d + neD + 2{n-'HD)ii + kdk)9D = k^s-ip, ( 4 . 8 ) 

in the conformal Newtonian gange. We follow the notation nsed in Ref. [21] with dots indicat¬ 
ing conformal time derivatives and 6d = (PD — P%)/p% and = V -dvD. Here, T~L = and 
%£) = Becanse the dark matter co-moving frame is not identical to the baryon one, 

additional bias terms proportional to (% — T-Ld) enter the above eqnations. The pertnrbed 
Poisson eqnation, in the baryon co-moving frame, is also modihed as 

5(j) -k 2>% {64> -k Td dTp) = —dvrGTv o? \5p + {Xd — 1) ^d\ ■ (4.9) 
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We implement these modifications in the CLASS code [24, 25] and show the effects of 
dark matter WEP breaking on the CMB power spectra in Fig. 7. For these plots, we fix 





Figure 7. The effect of varying rjjj = \d — 1 on the TT, EE, and TE power spectra for a transition 
redshift of zt = 10^ — 1. 


the transition redshift to be zt = 10^ — 1. For even a small deviation of Ad from one (or 
a tiny r/D = Ad — 1), dark matter WEP breaking has visible effects on the CMB power 
spectra. From Fig. 7, one can see that the main effect is to shift the location of peaks and 
troughs without changing the height of the peaks. This effect has some similarity to the effects 
from varying riA) except that the latter also changes the peak heights due to the late-time 
integrated Sachs-Wolfe (ISW) effect [37]. 


4.2 Constraints on Ad and zt 

Fig. 8 shows the constraint on the parameter ?7 d = Ad — 1- We see that the constraint on 
r/D becomes much stronger when the fluids decouple at high redshift (z.e., an earlier time). 
This is expected because the Friedmann equation is very sensitive to changes in the initial 
conditions. At high redshift, if Ad is not extremely close to one, an extraordinary amount of 
fine tuning is required to produce the observed universe. In the right panel, where we have 
fixed Zt = 10^ — 1, we show the two-dimensional contour in ryD and the current expansion rate 
Hq. The Planck three year measurement of Hq = (67.3 ± 1.3) kms“^ Mpc“^ [27], has some 
tension with the HST measurement of Hq = (74.8 ± 3.1) kms“^ Mpc~^ [28]. As one can see, 
by allowing a non-zero value of ryD, one can resolve the tension between these measurements 
in this dark matter WEP breaking framework. 
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Transition Redshift zt Expansion Rate Today [km s^' Mpc '] 

Figure 8. Left panel: 68% and 95% likelihood contours in the rjo = — 1 and zt plane. Right 

panel: 68% and 95% likelihood contours in the ijo = Ad — 1 and Hq plane. Here zt = 10^ — 1 is fixed to 
be a large value to provide a stringent constraint. Both plots use the dataset of Planck+Lensing+HST. 


5 Discussion and Conclusions 

In this paper we have updated the constraint on Newton’s gravitational constant G for all 
matter using the latest available cosmological data. We found a tension with the CODATA 
standard value for the gravitational constant which is independent of the specific combination 
of experimental data used in our analysis. The most stringent cosmological constraint on G 
comes from the combined dataset of Planck+ACT/SPT+Lensing+BAO+HST+BBN, and 
this result has a tension with the CODATA value at a significance of 2.2a. More insight on 
this discrepancy will come soon with the release of the 2015 Planck data [ 8 ]. If the best fit 
for G still tends high, this may be a hint for a new physics in cosmology and the need to go 
beyond the ACDM standard cosmological model. 

We also considered a traditional fifth force model to mimic the effect of equivalence 
principle breaking in the dark matter sector. In this type of model, the interaction between 
two dark matter particles is modified by the presence of an ultra light scalar, which mediates 
a long range attractive force. The main observable effect of this fifth force is to slow the decay 
of the gravitational potential at late times due to the expansion of the universe, adding power 
at low-.^ through the Integrated Sachs-Wolfe effect. Using the 2013 Planck data to constrain 
this effect, we found that the coupling strength of this fifth force relative to the gravitational 
force must be less than 10““^. This constraint is independent of the mediator mass for the 
region of parameter space which does not modify the cosmological background evolution. 

Finally, we considered true WEP principle breaking in the dark matter sector (be., an 
explicit difference between dark matter gravitational and inertial mass). General relativity, 
being a geometrical theory, encodes this equivalence by construction. As such, to break 
the WEP in the dark matter sector we follow a phenomenological approach, motivating the 
modifications to the cosmological equations from post-Newtonian and non-relativistic fluid 
dynamics where this principle can be easily broken. Providing a fully consistent theory 
that encodes this phenomenology is beyond the scope of this work. The main thing to be 
learned here is that such modifications drastically modify the evolution of the cosmological 
background. The observable effect, even for extremely small degrees of WEP breaking, is to 
cause a visible, dark energy like shift the location of the peaks in the CMB power spectrum. 
Within the framework of this phenomenological model, we found that the ratio of the dark 
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matter gravitational to inertial mass is constrained to be unity at the 10 ® level when the 
two fluids decouple at times earlier than approximately the recombination time. 
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A Linear Perturbation Equations for the Two-Fluid Description 

For our two-fluid description, we have two independent background evolution functions: a{t) 
for ordinary baryons and auit) for dark matter. We use a{t) as a clock such that the linear 
perturbation for ordinary baryons will be unchanged, except that the common gravitational 
potential is also sourced by dark matter. 

Following the Newtonian approach in Mukhanov [17], the basic equations that govern 
the baryon and dark matter perturbations are: 

+ Va;(/Jv) = 0, (A.l) 

d'v 

+ (v • Va:)v -|- x<i> = 0 , (^-2) 

Vl(t> = 4.TTGNP, (A.3) 

where are continuity, Euler and Poisson equations, respectively. For the dark matter held, 
one can separate the helds into the zeroth and leading order as 

PD = Pd + ^Pd , = ^D + , (j) = (l)° + 6(j) . (A.4) 

Substituting the perturbations into the above three equations, we obtain the following equa¬ 
tions 


—^-h + Va;((5/?D • Vj)) — 0 , (A.5) 

+ i^D • ^x)S^rD + • v.)v^ + VxScI) = 0 , (A.6) 

'V'^64> = AttGn [Spb + Xd ^Pd] , (A-7) 

The last Poisson equation should be modihed when including radiation and dark energy 
densities. 

For the continuity equation and using the background held x, we have 3H o = 

Vx • v^). Here, the Hubble parameter is dehned as Hd = The zeroth order continuity 

equation has the canonical relation dp^/dt = —3Hjjp^. Dehning the dimensionless quantity, 
6d = dpn/p^, we have the hrst-order continuity equation as 

^ - (V^ • 6d + Vx- + VxidD • V?,) = 0 . (A.8) 
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In order to only use the ordinary baryon scale factor as a measure for time, we transform 
the dark matter equations to the baryon co-moving frame with x = a(t)q. The relation 
between the time-derivative’s in the universe and co-moving frames is 


d 


d 


dt 


={m) -K-V.) = 


= fl'] - (JL 


dt 


\Hd 




V, = -Vg = -V,(A.9) 
a a 


where we have used the relations x and x to have Hd. So, in 

the baryon co-moving frame, we have 


^ + -V • hvB = 0, (A.IO) 

ot a \ Hd J a 

+ - f 1 - {vd • V)(5vd -h Hd 6^rD + -V(5(/) = 0, (A.ll) 

at a \ Ho J a 

V^(5</> = 4:TtGn c? [/3j (5f, -|- \o p% (^d] • (A. 12) 


Dehning the divergence of the velocity held as 6o = Vhv^) and changing the ordinary-time 


derivative to confermal-time derivative, we have 

+ (l - + 9d = 0 , (A.13) 

^D+{l- V • [(v?) • V),5 vz5] + HdOd + = 0, (A.14) 

= AttGn [pI) hfc -|- Xo p% So] ■ (A. 15) 

Here, we have used H. = aH and Ho = olHo- After some algebra, one has 

V • [{v% • V)hvB] = SnoBo + v?) • V^Svo + {.A • V)0Z) (A. 16) 

= 3no9o + 2iv%-V)9o. (A.17) 

In the second line of the above equation, we have used V^A = V(V • A) — V x (V x A) and 
assumed curl-less velocity perturbations (he. V x = 0). We now dehne the operator D 

^=(l-^)(^?5-V). (A.18) 

The continuity and Euler equations become 


So + i^So + 9o = 0 , (A.19) 

9o + {4:no-3n + 2'D)eo + X7^S^ = 0. (A.20) 

When Ho = ^ (*-e. normal cosmology), the operator T> is exactly zero and the equations 
reduce to the equation in the standard cosmology [21], The appearance of this new operator 
indicates the breaking of equivalence principle. Since we chose the frame to be the baryon co¬ 
moving frame, the operator T) appears as an artifact in the dark matter equations to convert 
from the dark matter co-moving frame to the baryon co-moving frame. Using the Fourier 
transformation in Appendix B, we have the corresponding equations in momentum space as 


So + {H — Ho) {3 + k dk) So = 3S^ — 9o , 

9o + H9o + 2{H — Ho) (1 + ^ dk) 9o = S'ljj, 
k"^ Scj) + 3H{S^ + HSijj) = —AttGn [p^ Sb -I- Xo p% h/j] • 
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(A.21) 

(A.22) 

(A.23) 




In the Newtonian gange and with the assumption of no shear, we have 54> = Si/j. In the right- 
hand side of the hrst eqnation of Eq. (A.21) and in the left-hand side of the third eqnation of 
Eq. (A.23), we have also added an extra term to acconnt for the geometry part from general 
relativity, which is beyond this post-Newtonian description. 

B Fourier Transformation of the New Operator V 

Let ns consider the Fonrier transform of this new operator acting on some fnnction /(x) = 

f{a ■ q) 


T{Vf) = I = ( 1 - ^ • V/). 


'V 


i-Al 




Td-D ) 


'V 


Integrating it by parts, one obtains 




[ e-^^-^{v%-n)fdS- [ A/V- 
JdV Jv 


(B.l) 


(B.2) 


where the snrface integral vanishes when integrating over all space as long as / decays fast 
at inhnity. So we have 

= (^ " ^) X • V?)) + V • v^] 

d^q • q) -h SUd] 


= f--l 

\nD 

= {n- Ud) 


IV 


3.A(/)+k-Vk / d^qfe 

Jv 

= {n- Ud) [3 + (k • Vk)] Hf) • 


—zkq 


(B.3) 


Where we have nsed = Hu-x. = %£)({. 
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